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We  give  an  explicit  formula  which  describes  the ‘solution  of  the 


problem  of  the  vibrating  string  with  a  plane  obstacle.  This  formula 


allows  us  to  prove  continuous  dependence  on  the  data;  a  regularity 


result  is  given.  We  prove  some  results  on  the  convergence  of  the 


penalized  problem,  and  give  a  numerical  scheme. 


A  few  results  are  given  without  the  requirement  that  the  obstacle 


be  plane. 
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SIGNIFICANCE  AND  EXPLANATION 


The  purpose  of  this  paper  is  to  continue  the  study  of  the  following 
problem:  a  string  is  constrained  to  remain  on  one  side  of  a  material  obstacle. 
We  consider  its  transverse  vibrations,  and  we  assume  that,  when  the  string  hits 
the  obstacle,  no  energy  is  lost.  It  was  proved  in  a  previous  paper  that  when 
the  obstacle  is  concave  i.e.  has  the  form  of  a  bowl  (with  the  plane  obstacle 
as  a  limiting  case) ,  then  there  exists  a  unique  solution  for  a  given  initial 
position  and  velocity  of  the  string.  We  prove  here  that  when  the  initial  posi¬ 
tion  and  velocity  are  slightly  changed,  the  solution  is  only  slightly  changed. 
Moreover,  if  we  replace  the  rigid  obstacle  by  an  elastic  one,  if  the  obstacle  is 
plane,  and  if  the  initial  data  are  sufficiently  regular,  then  the  solution  we 
obtain  is  close  to  the  solution  of  the  original  problem.  We  also  give  a 
numerical  scheme  for  computing  the  solution. 

This  subject  is  a  first  approach  to  considering  problems  of  mechanical 
vibrations  with  unilateral  constraints;  a  number  of  elliptic  and  parabolic 
problems  with  unilateral  constraints  have  been  solved,  but  very  few  hyperbolic 
problems,  i.e.  problems  describing  vibrating  mode  s  such  as  the  above  obstacle 
problem,  have  been  tackled  successfully. 
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NEW  RESULTS  ON  THE  VIBRATING  STRING 
WITH  A  CONTINUOUS  OBSTACLE 


A.  Bamberger^ 

M.  Schatzman 

I .  INTRODUCTION 

1.1.  Presentation  of  the  problem  and  the  results. 

This  paper  aims  to  give  some  new  results  on  vibrating  strings  with  obstacles.  The 
model  is  the  same  as  in  [1] ,  but  as  it  appears  necessary  to  elucidate  several  points  of 
the  modelization  which  was  exposed  there,  we  shall  give  it  from  the  beginning. 

We  consider  the  small  transverse  vibrations  of  a  string  that  is  constrained  to  be 
on  one  side  of  a  material  obstacle.  Let  the  transverse  displacement  at  time  t  of  the 
material  point  of  the  string  with  coordinate  x  be  denoted  by  u(x,t).  If  the  string 
were  free,  i.e.  if  there  were  no  obstacle,  then  u  would  satisfy  the  wave  equation 

□  u  =  u  t  -  u  =  0  . 

tt  XX 

We  assume  that  the  obstacle  has  position  ^(x)  .  We  translate  the  requirement  that  the 
string  stay  on  one  side  of  the  obstacle  into  the  inequality 

(1)  u(x,t)  >_  *  (x)  V  x,t  . 

When  the  string  does  not  touch  the  obstacle,  its  motion  satisfies  the  wave  equation, 
and  thus 


(2)  suppDu  c  {  (x,t)/u(x,t)  =^t(x)}. 

We  require  that  the  string  does  not  stick  to  the  obstacle;  this  can  be  translated  as 

(3)  Qu  >_  0  , 

which  means  that  the  obstacle  does  not  exert  a  downwards  force  on  the  string. 

Notice  that  (3)  is  essentially  equivalent  to  subsonic  propagation  of  interactions. 
To  see  this,  let  t  =  c  (x)  a  curve  which  separates  a  region  *<  of  the  half  plane 
3R  *  (0,-)  in  two  open  regions  *f+  and  where  Hu  vanishes.  Suppose  that  - 
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and  u 


are  sufficiently  smooth,  and  that 


(4) 

(5) 


u‘  (x,o  (x) )  =  v'  (x) 

u“(x,t)  >  ^(x)  v  (x,t)  Hi. 


Then  wc  can  compute  Du  in  the  sense  of  distributions,  with  ■,  a  test  function: 


(8) 


<  l  J  u ,  V  '  =  -  <  - -  t  “T  /  +  '  T. —  ,  ~~z~  >  = 


/[[  w-  j(x'0(x))  + 


uT  ) 

•  X  I 


(x,  (x)  )  *  (x)  I  .  (x  ,  ( X )  )  dx. 


I 


Relation  (4)  can  be  differentiated,  with  respect  to  x  ,  and  implies 


(7) 


3u 


\ 


- - -  (x,0  (x)  )  =  -o'  (x)  —7 - 7—  (x,  (x) )  . 

ax  ax  |  |  .it  *t  I 

Introducing  (7)  into  (8),  we  get: 


<  □  U  ,  t;  >  -  f  |  j  (x,<:  (x) )  (1  -  •  ,2  (x) )  .  (x,  :  (x)  )dx. 

But  hypothesis  (4)  and  (5)  ensure  that 

(x,o(x))  _  0  and  (x,.(x))  _0  . 

Therefore,  Du  is  non-negative  if  and  only  if  | «> •  j  is  almost  everywhere  smaller  than  1. 

It  is  not  enough  to  suppose  that  conditions  (1),  (2)  and  (3)  arc  satisfied,  as 
nothing  has  been  said  of  the  evolution  of  the  energy  of  the  string  during  the  collision 
with  the  obstacle. 

The  hypothesis  that  will  be  made  is  that  the  energy  is  conserved.  This  requirement 
should  be  analysed  from  a  mathematical  point  of  view  as  follows:  the  condition  must  be 
local,  because  the  propagation  properties  of  hyperbolic  equations  suggest  it,  and  it 
must  be  satisfied  wherever  in  the  x,t  half-plane  the  free  wave  equation  is  satisfied. 
Thus,  multiplying  by  — the  relation 
(8)  I  '■  u  =  0  on  *f, 

where  rtf  is  an  open  region  such  that  (8)  in  satisfied,  we  obtain  a  rc  lation  in  divergence 
form 


Thu  operations  by  which  wo  deduce 
square- i ntey ruble  in  U  •  (<),). 


j  fn  _ai  *u  '• 

* x  i  “  <t  x  ; 

C3)  out  of  (H)  are 


in 

va  1  i  d 


if 


•u 

T 


and 


art  locally 


we  shall  prove  in  paragraph  II  the  following  explicit  formula  in  the  case  of  the  plane 


obstacle. 

Let  w  be  the  free  solution  of  the  wave  equation 
□  w  =  0 

w(x,0)  =  uQ(x)  , 
wt(x,0)  =  u^tx). 

r.et  the  obstacle  be  =  0  ,  and  let  the  backward  wave  cone  be 

T-  d§f  {<x',f)/0  <_  f  <_  t  -  (x-x*  J  >. 
x  ft 

Let  us  denote  by  r  the  negative  part  of  a  number 
r  =  sup(-r,0)  . 

Then  the  solution  of  the  problem  (1)  -  (3)  and  (11)  is  given  by 

u(x,t)  =  w(x,t)  +  2  sup{  (w(x'  ,t' ) )  /(x’.t')  e  T~  }. 

x,t 

This  formula  shortens  considerably  a  previous  proof  [6]  of  continuous  dependence  on 
data,  and  is  the  key  for  the  numerical  scheme  studied  in  paragraph  III.  We  shall  give 
ir  paragraph  IV  a  regularity  theorem  in  spaces  of  bounded  variation,  in  the  case  of  a 
general  concave  obstacle. 

In  paragrapn  V,  we  shall  consider  the  functions  u^  which  solve  the  problem 

Du)  -  X  <UA  -  *>'  ■  0 
u^ (x ,0)  =  uQ(x) 

(x.0)  =  Uj^tx) 

In  the  first  part  of  this  paragraph,  we  shall  prove  a  weak  convergence  result,  which  doe 
not  depend  on  the  shape  of  or  on  the  regularity  of  the  initial  data.  The  limit  func 

tion  will  satisfy  a  set  of  energy  inequalities  instead  of  (11). 

duo 

In  the  second  part,  we  shall  assume  that  the  obstacle  is  plane,  and  that  -njj-  and 

Uj  are  locally  of  bounded  variation.  Then  the  solution  of  (16)  converges  strongly  in 

H,  (IP  ■  IP*  )  ,  and  its  limit  is  the  unique  solution  of  (l)-(3>  and  (11). 

loc  ^ 
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1.2.  Notations  and  summary  of  previous  results. 

We  shall  use  throughout  this  paper  the  following  notations  and  definitions: 
V  is  the  set  of  functions  u  such  that 


(17) 


/a  ( |u  (x,t)  | 2  +  |ufc  (x,r)  |2)  dx  <_  C(a,b)  <  +«■  Va,b  v  t  <  b 


w  is  the  free  solution  of  the  wave  equation: 
□w  =  0  , 


(18) 


w (x,0)  =  uq(x) 
wt(x,0)  =  u^ (x)  . 


E  is  the  set 


E  =  {  (x,t)/w(x,t)  <  (x)  } 

I  is  the  domain  of  influence  defined  by 

(19)  I  =  u{T+  /(x,t)  £  E) 

where  T+  is  the  forward  wave  cone  {(x'.rU/t1  >  t  +  |x-x’|},  and  its  boundary  is 
x#r 

given  by 

31  =  {  (x,r)/t  =  t (x) } 

where  i  is  Lipschitz  continuous  with  Lipschitz  constant  1  . 

(see  [1],  Proposition  II. 3  for  the  proof  of  this  claim). 

The  backward  wave  cone  T  is  {(x',t')/0  <_  t'  <_  t  -  |x-x'  |  }. 

The  characteristic  coordinates  £  and  r\  are  given  by 

x+r 


(20) 


(21) 


_  x+r 

5  =  7T 


n  = 


with  the  notation  z(£,n)  =  z  (-^=^-  ,  ^  for  aH  functions  of  two  variables  x  and  t. 

i  -> 

We  shall  call  problem  (P^)  the  following  problem:  given  uQ 
satisfying  the  compatibility  condition 


H.  UR),  u, 
loc  1 


L,  (JR) 
loc 


(22) 

(1) 

(2) 

(3) 


u  (x)  >_  (x) 

n^(x)  >_  0  a.e.  on  (x|u0(x)  =  ^  (x) } 
find  u  in  V  such  that 

u  >  v5 


supp  Du  c  { (u^tj/utx.t)  =  *£  (x)  ) 
O  u  >  0 
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i  Z  2  'j 

—  (u  +  u  )  +  -r—  (-2u  u  )  =0  in  the  sense 

Uo)  1  x  c  Jx 

i  of  distributions  in  E3  x  ir+ 

1  u  (x,0)  =  u  (x) 

(23)  j  .  0 

'  -£•  (x,0)  =  (x) . 

The  precise  statement  of  the  results  of  existence  and  uniqueness  in  [ 1 ]  is  as  follows: 
Theorem  0.  The  Problem  P  possesses  a  unique  solution  u  if  is  non-negative. 

Moreover,  this  solution  u  is  the  unique  solution  of  the  linear  problem 


| { (x,t)/t  f  t (x)  } 

—  (x,  t  (x)  +  0)  =  -  (x,t(x)  -  0)  a.w.  on  jx/r  (x)  •  0  s.  |  I  1  (x)  ;  -  1 : . 

-jt  ot 

u (x, 0)  =  uQ(x) 

(x,0)  =  u1 (x ) 


If  u  is  the  measure  defined  by 


<  u,  v  >  =  -2  /w  (x,  t  (x) )  (1  -  i  '  (x)  2)  '(j  (x,  t  (x)  )dx  , 


then  the  solution  of  (19)  is  given  by 


u  =  w  + 


where  &  is  the  elementary  solution  of  the  wave  equation  with  support  in  the  positive 
light  cone: 

I  fe  =  'j  on  {  (x,t)/t  |x|  > 

(27)  j 

^  &  =  0  elsewhere. 

It  will  be  useful  to  consider  the  problems  P  .  ,  which  are  -just  P  restricted  to  the 

x,t  J 

backward  wave  cone  T  ,  with  initial  data  given  on  (x  -  t,  x  +  t).  Clearly,  u  is 
x ,  t 

a  solution  of  P  if  and  only  if  it  is  a  solution  of  P  for  all  x  IP.  ,  t  ">  0  . 

x ,  t 

The  first  result  on  the  convergence  of  the  penalty  method  for  the  string  with  an  obstacle 
was  proved  by  A.  Bamberger  [21. 

An  explicit  formula  for  the  string  with  a  point  obstacle  was  obtained  by  L.  Amerio 
in  [3]  and  by  M.  Schatzman  in  (4],  with  a  different  argument. 
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IX.  the  explicit  formula,  continuous  dependence  on  the  initial  data 


II. 1.  The  explicit  formula  for  the  infinite  string. 

In  the  case  of  the  zero  obstacle  (and  more  generally,  the  plane  obstacle),  the  solu¬ 
tion  of  (P^)  can  be  expressed  by  an  explicit  formula.  We  denote  by  r  =  sup (r,0)  the 
negative  part  of  a  number. 

Theorem  1 .  The  unique  solution  of  P  when  #  =  0  is  given  by 

(28)  u(x,t)  =  w(x,t)  +  2  sup  (w(x',t')] 

(x*,t')€T' 

X  r  u 

Remark  2 .  If  the  obstacle  is  plane,  i.e.,  if  sf(x)  =  ax  +  S,  then  (28)  can  be 
generalized: 


(29) 


u(x,t)  =  w(x,t)  +  2  sup  Iwfx’.t1)  -  V’(x’)  ]  . 


(x*  ,f)€T 


x,t 


To  deduce  (29)  from  (28)  it  is  enough  to  consider  u  -  ,  and  notice  that  C  C  =  0.  ■ 

The  proof  of  Theorem  1  comes  in  several  steps. 

The  first  step  is  the  following  result: 

Lemma  3.  The  set  where  sup{  [w(x*  ,t* )  ]  /(x’.t1)  e  T..  does  not  vanish  is  the  interior 

x ,  t 

of  the  domain  of  influence  I. 

Proof.  If  w(x',t’)  <  0  for  some  (x',t')  in  the  backward  cone  T„  _  ,  then  (x,t) 

x ,  t 

belongs  to  the  forward  cone  T+,  ,  the  vertex  of  which  is  in  the  interior  of  r  .  Thus 

x  ,t 

(x,t)  is  in  the  interior  of  I  .  Conversely,  if  (x,t)  belongs  to  the  interior  of  I  , 
then  there  exists  a  point  (x',t*)  in  the  interior  of  E  such  that  (x,t)  belongs  to 

the  interior  of  T+.  We  can  choose  this  (x'.tM  such  that  w(x't’)  be  strictlv 

x'  ,t ' 

negative,  because  the  set  of  (x',t‘)  such  that  w(x',t’)  <  0  is  dense  in  the  interior 

of  E  .  Therefore,  sup  [w(x',t')]  >  0  .  ■ 

(x'  ,t')rT_  . 

x,t 


(30) 

Then 


(31) 


Let  us  define 

k(x,t)  = 
,  thanks  to  Lemma  3 , 
u(x,t)  = 
u(x,t)  = 


inf  {w(x 
we  have , 
w(x,t) 
w(x,t)  - 


,f)/(x',t<)  ,  T"  }  . 

xtz 

if  u  is  defined  by  (28) 
for  t  <_  r  (x) , 

2k(x,t)  for  t  ^  i(x)  . 
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Let  w =  f(;')  +  g(n),  where  f  and  g  are  in  (IR)  .  From  (35),  we  deduce 

loc 

f('")  +  g('i')  '  0  for  )  and  from  (34),  f(  r)  +  g(n  )  >0  for  ' .  r  ■  •' !  . 

l—  l—  “i  i—  l—  “  l 

As  we  must  have  fC.)  +  g('i.‘)  =  0  =  f  ( "  • )  +  g  ( n ! )  =  0  ,  by  definition  of  U,  ,  ■'.!  , 

11  11  1  1  1  1 

and  ■.!,  then 

i 

(37)  f  (.'.>=  f(-!)  -  f(i)  V  ■ 

i  i  —  l  i 

Similarly , 


(38) 


On  C 


(39) 


g(:ij  =  g(nM  _g(n)  V  n  [nl.n.J  . 

the  complement  of  the  set  u  [  -  ,  f],  we  have 

l  i 

l 

f 1  ( 5 )  '  0  a.e.  . 


The  simplest  way  to  see  this  is  to  notice  that  Y  is  one  valued  on  C  ,  and  that 
HZ)  +  g  (Y  (i)  )  =  0  on  C 
f  ( C  * )  +  9  (Y  (i) )  0  for 

Let  us  evaluate  now  f(;)  +  g(-;).  Suppose  first  that  X ( n)  =  Y  On)  is  one  valued. 
Then 


x(n) 


(40) 


f(f)  +  g  (n)  =  g  ( r  i )  +  f  (x  (n) )  -  l 


lc  f(-*)df  - 


(f  (min  (tpx(n) ) )  -  f  (max  (>:,£))  ) 


where  the  summation  is  extended  to  the  indices  such  that 
We  have : 


intersects 

i  l 


(i,X(h)l. 


g  (n)  +  f  (x  (n) )  ^  o  . 
x(n) 

j  1  f  •  < ' * ) d - •  ^  0  by  (39) . 

As  X  ( ■ )  is  one  valued,  it  is  not  contained  in  the  interior  of  an  interval  {  ■'  ."'l. 

i  i 

Thus 


min  if  [•".  , r'  1 

i  i  li 

the  corresponding  term  in  the  sum  vanishes.  For 

f  (••.!)  -  fC)  , 

i 


t ' .  X  ( •  )1  f  $  and,  it  ■'  \  (y.-'M, 

in  the  term  in  the  sum  is 

l  i 


which  is  net  positive,  by  (3.7).  Therefore, 
X(-i).  If  we  suppose  that  X(-  )  =  (  '. 


3  3 


y  (  ,)  +  f  (:!  -  /  ]  lpf  •  <  -I  * )  tl '  • 
and  the  result  still  holds,  i  .  o  .  : 


the  expression  (40)  is  non-negative 

] ,  we  have  to  study  the  expression 

( f  (min  ( '  !  ,  ' )  -  f  {max  ('.,*)) ] 
i  )  i 


for 
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(41)  w(x,t)  2.  0  for  t  <_  T  (x)  . 

Thanks  to  (41)  ,  we  may  redefine  k  as 

k(x,t)  =  inf  {w  (x1  ,  t '  )/t '  <_  t  -  |x-x'|}, 
or  still,  in  characteristic  coordinates, 

(42)  kU,n)  =  inf{f(£’)  +  g(n')/C*  <  (  M'  i  n). 

Then,  it  is  immediate  that 

(43)  k(C,n)  =  inf  { f  (£')/£ 1  <_  O  +  inf{f(n')/>r  1  n), 

which  proves  the  claim  of  Lemma  4.  ■ 

We  shall  now  prove  that  u  ,  defined  by  (28)  satisfies  the  transmission  condition 
across  the  line  of  influence. 

Lemma  5.  If  u  is  defined  by  (28),  then  almost  everywhere  on  {x/|t*(x)|  <  1}, 


(44) 


1)11  (x,t(x)  +0)  =  -  ~  (x,t(x)  -  0). 


3t 


Proof.  Let  A  =  { x/ | T ' (x) |  <  1}.  Then,  almost  everywhere  on  A  ,  by  Corollary  A. 2  of  [1) 
(45)  wx<x,t(x))  and  wt(x,i(x))  exist. 

Let  x  be  a  point  satisfying  (45) ,  and  let  us  denote 

w^  (x,x  (x) )  =  a  ,  wfc(x,T(x))  =  b  ,  T  *  (x)  =  m  . 

Then, 

a  +  mb  =  0  ,  b  £  0 

and 

w  (x 1  ,  t 1 )  =  a(x'-x)  +  b(t*  -  t(x))  +  t(x'-x,  t'  -  t  (x) ) 


where 


lim  e(r  x) 

|  r  |  +  |  s  |  ->0  FTR=  ° 


We  have 

inf{a(x'-x)  +  b(t' 
and  therefore , 


T (x) )/(x* )  <  T  }  =  b (t  -  T (x) ) 

X  f  t 


! 


b  (t 


T  (x)  )  -  supi  I  C  (x'-x,t’-t)  j/l  <X'  )  ^  t'  <_  t 
k(x,t)  <  b(t  -  r(x))  +  |f(0,t  -  t  (x) )  |  . 


x-x' 


(46) 


■  UUW:LIIW«F, 


As 


I  t  '  (x) I  <  I  ,  we  have 


lira  [sup{ | e (x'-x,t'-t)/T (x1 )  <_  t*  <_  t  -  |x-x'|/t  -  T  (x)  ]  =  0 
t+T (x) 

and  we  deduce  from  (46)  that 

u(x,t)  -  u(x,t(x))  ,  ,  .. 

lira  - T  — r-j - -w(x,t(x>) 

t  *  T  (x)  W 

under  the  assumption  (45) .  • 

Conclusion  of  the  proof  of  Theorem  1. 

Lemmas  3,  4  and  5  imply  that  the  function  u  defined  by  (28)  solves  the  linear 
problem  (24) ,  up  to  the  condition  u  e  V  .  Therefore  it  remains  to  check  this  last  condi¬ 
tion.  If  we  take  into  account  the  formula  (43)  ,  let  us  show  that  k  is  in  V  . 

We  know  that  f  is  in  H?’  (3R);let 

loc 

f (5)  =  inf{f (?•)/£•  <  C). 

Then,  we  can  compute  the  derivative  of  f(£)  almost  everywhere: 

f’(0  -  0  if  f(?)  >  f(C)  or  if  f'U)  >  0 


(47) 


f'(C)  -  £•(£)  if  f(?)  =  f(e)  and  if  f’(S)  <  0 


We  deduce  from  (47)  that  f  is  in  H*  (IR)  .  Similarly,  g  is  in  H?’  Cm).  The 

loc  loc 

function  k  which  can  be  written  as 


k(x,t)  =  f(^)  +  g( 


V2  • 

will  therefore  be  in  V  ,  i.e. 

3  ,2 


~7T 


/  (|kx(x,t)  |  +  |kt(x,t)  |  )dx  <_  C(a,b)  V  a,b,  V  t  <_ 


and  thus  u  is  in  V  . 
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L,  <  IF.)  satisfying 
loc 


IX, 2 .  Continuous  dependence  on  the  data 
Corollary  6.  The  map  (Ug,^)  -»  u  which  to  an  element  of  H^^dR) 

the  compatibility  condition  (22)  associates  the  solution  of  is  continuous  from 

1  2 

H  (R)  x  L,  (IR)  equipped  with  the  strong  topology  to 
loc  loc 

W*'P([0,+“);  L?  (IR))  n  I*  ([0,+“>j  H*  (IR)) 
loc  loc  loc  loc 

equipped  with  the  strong  topology,  for  all  finite  p  . 


Proof.  We  have  at  once  the  continuity  from  (IR)  *  (IR) 


to  C  (®  <  »  )  . 

The  topology  on  W^P  ( [0,+«°) ;  L^dR)!  n  LPnf,(  [0,+°°)  OR) )  is  defined  by  the 


semi-norms  for  A,B  >  0 


loc 


q  (u)  =  | u (0 ,0 ) |  +  (/B  </A ( u2  +  u^) (x,t)dx)P/2)1/P  . 

ABP  j  0  -A  X 

The  topoloqy  on  H  (IR)  x  L,  (IR)  is  defined  by  the  semi-norms  for  A  >  0 
loc  loC  .  ,du  ,2 


PA<U0'U1)  =  |U0<0)|  +  ( 

-A 


_ 0 

dx 


+  lu, |2)dx)1/2 


It  has  been  proved  in  (1) ,  paragraph  IV. 2  that  for  solutions  of  P  with  zero  obstacle 


(48) 


Therefore 


iff  tx.t) I  =  iff  <*+t,0) | 

I  3u 


||£  (x,t)|  =  (x-t,0) 


(49) 


Let  q 


r  duvr  + 1 u  i  > (x,t)dx  =  r  <ki  + 1«  r> <x,t>dx  < 

-a  4  n 

-  ,  ,du  1 2  ,] 

-  >  -ST  +  IU1 1  d*  v  >  0  . 


AB« 


A+t 

/ 

-A-t 

be  the  semi-norm 


qAB»(V> 


r(0,0)  +  ess  sup  (  /A(|v  |  +  | v  |  ) (x,t)dx) 

t<  [0,B]  -A  * 


1/2 


Then,  (49)  implies 

<50>  qAB~(U>  -  PA+B(U0  +  V’ 

If  (Ug,uP)  is  a  sequence  of  initial  data  satisfying  the  compatibility  condition  (22) 

and  converging  to  (u „,u  )  in  H?'  (IR)  x  (IR)  ,  then,  as  a  consequence  of  (52), 

0  1  loc  loc 


(51)  un  -  u  in 

and  moreover,  (48)  implies  that 


H^oc(IR  x  ir+)  weakly 
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(52) 


n 


(x,t)  “dx 


V  t ,  A  ■  0 


—  (x.t)  -dx  .  (  - 

—A  -A 

a  'un  ,  .  2  fA  •! 

j  — —  (x,  C )  ,  dx  •  1  — 

-A  -A 

3atnering  (51),  (52)  and  (53),  we  obtain 


(53) 


(5-1) 


u 


r (x , t )  dx 

-  (x.t)  “dx 

1  , 
h  ,  ( n- 

loc 


V  t  ,s 


B  ) 


Thanks  to  rubini's  theorem,  one  lias  from  (54) 


st  rongly . 


(55) 


(u  ( - ,  t )  ,  u  ( - , t ) )  •  (u  (*,t),u  (  *  , t )  ) 

X  t  X  t 


in  (I'  (BO)' 

loc 


strongly,  for  almost  all 
The  relation  (55)  together  with  the  estimate 


'*AB-,U  > 


•SUI  A»B(U 
n 


n\ 

ui) 


imply  that 


H.  (IP) ) 
loc 


converges  to  u  in  the  space 


Remark  7. 

T'ne  mapping 

(uo,ui>  ’  u 

i: 

s  not  continuous 

Lloc(I°'  + 

■ )  ,-hJ'  (IP)  ) 
loc 

which  is  the 

space*  V  defined 

Take 

for  instance 

the  sequence 

of 

initial  data 

n  i 

uo  =  1 

n  n+ 1 

-  3i =  ~  • 

As  these 

do  not  depend 

on  x  ,  the  solution  of  p  i’ 

n  .  . 

u  (x,t)  = 

1  i  -  n±i  t 

{ 

2iit-l 

\  n 

if  t  -~9 

n+ 1 

if  t  -2—, 

—  n+1 

with  the 

limit 

1  -  t 

if  t  1 

u(x,t)  = 

1 

t  -  1 

t 

if  t  _  1  . 

Then  we  may  calculate 

qAB  _  U) 

,  11 

qAB-  <U 

-  u)  =  0 

if 

*  ^ 

/  n 

qAB<U 

-  u)  -  2.~ 

if 

*  • 

Thus  if 

B  1  *  qAB 

(uP  -  u)  doe 

not  tend  to  zi  r 

<}'  r  (('),+  •)  ;lT  (  D-  )  ) 

loc  loc 


loc 


loc 


i.;  (  b-  ) ) 

lot- 


inf  ini  tv. 
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III. 3.  Application  of  the  explicit  formula  to  the  finite  string  with  fixed  ends. 

The  explicit  formula  (29)  will  allow  us  to  give  a  simple  construction  of  the  solution 
of  P^  where  P^  is  the  problem  of  the  vibrating  string  with  fixed  ends,  and  obstacle 
v  -  -K  <>  .  The  only  modification  with  respect  to  P  ,  is  that  we  shall  require  u  to 

be  in  the  space  L  (u  ,T;  hj  (<» ,  L) )  W*'  (u  ,T;  \?  (U ,  L) )  for  all  TO. 

In  fact  u  will  be  in  the  space 


L 11^(0, L>)  W1,'((U,')i  L2(0,L)) 

because  we  can  integrate  (10)  on  any  rectangle  (0,L)  *  (0,T),  and  we  get  the  energy 
equality  for  arbitrary  times  T  : 


(56) 


/L(  ,Ux<x,T)  ,  J  +  |ut  (x,T)  j  2 )  <Jx 


. 

i  dx 


Let  us  define 


e 


Then 


f  ( 

!/Ll 

"0  ■ 


^0 

dx 


1  1 1/2 


1 


dx 


,u(x,t)  -  u (o ,t)  ; 


and  similarly 


| f*  u  (x‘  ,t)  dx'  [  e>x 
0  X 


ju(x,t)  -  u(L,t) ;  _e>L-x  . 


Let  .  =  ,  Then 

e‘ 

(57)  V  t  ('),•),  V  x  [(),,)  !L  -  .,L1  ,  u(x,t)  -K 

and  l.u  cannot  be  supported  in  the  strips  ([0,.)  (1.  -  .  ,Ll)  •  f  0 ,  • )  . 

the  initial  conditions  u(),  u^  to  the  interval  l-t,L  +  ,]  by: 

u  (— x )  -  -u.  (x)  if  x  l  =  0,1. 

li 

uMLtx)  =  -upi.-x)  if  x  10,  I  ,  i  =  0,1. 

Then  tlie  corresponding  fret  solution,  w  is  defined  on  the  cone  ^ 

~2  1  3" 


property  that 


Let  us  extend 


,  with  the 


w(x,t)  *  o  <)  _  t 
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Let  u  be  defined  on  T 


by  (29);  then  for  x  =  0,  t  <  a 


2  '  2  +  a 

u(0,t)  *  w(0,t)  +  2  sup  {[w(x',t')  +  KJ  }. 


But,  w  (0,t)  =  0  ,  and  T  is  included  in  the  strip  I-a,aJ  *  10,“),  so  that  w 

o  ,t 

on  this  strip,  and  thus  u(0,t)  =0  on  [0,a].  Analogously,  u(L,t)  =0  on  [0,  .] 


Therefore,  (5.8)  defined  the  solution  of  P,  on  T 

f  L  L 


n([0,L)  «  [0,”)), 


2  '  2  +  * 

Let  us  define  by  induction  the  solution  of  P.  on  the  region  R  given  by: 

f  n 

f  R„  =  {  ( x , t )  •  (0 ,L]  x  [0,“)/  \  +  (2n-l)i  -  |x  -  <_  t  <_ 


<_  j  +  (2n+l)a  -  |x  -  j| )  ). 


We  shall  denote  by  a  the  function 
'  n 

f  un(x)  =  j  +  (2n-l )  a  -  jx  -  j  |  if  x  .  [0,L! 

I  o  (x)  »  o  (— x)  o  (L  +  x)  =  o  (L  -  x)  if  x  ■  [ 0 ,  .  ]  . 

n  n  n  n 

Suppose  we  know  u(x,on<x))  for  X  •;  [0,L],  Let 
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w  (x,o  (x) )  *  u(x,o  (x) ) 
n  n  n 

if 

X  • 

10, L] 

w  (x,o  (x))  =  -u(x,o  (x> ) 
n  n  n 

if 

X  • 

[0,  I 

w  (x,o  (x) )  =  -u(2L  -  x,  o  (2L  -  x) 
n  n  n 

if 

X  V 

IL;  L  + 

^  □  wn  =  0 

The  function  wn  is  defined  in  the  region 

°n  (x>  i  i  (2n  +  “  +  J  *  |x  -  j  | , 


x  <  L  +  i. 


Let  us  notice  that  the  sytmetry  of  the  initial  conditions  in  (60)  implies 

(61)  w  (0,t)  =  w  (L,t)  =  0  for  (2n-l)j  <  t  <  (2n+l)i. 

n  n  —  — 

Moreover,  as  u  (x,a^(x))  2.  "K  <  V  x  f  {0,L]  ,  and  as  u  satisfies  the  energy  condition 
(10) ,  we  shall  have 

w  (x,t)  >.  -K  for  o  (x)  <  t  <_  (2n+l)or  +  -  |x|  or  for 

(62)  •  "  n 

°n*x'  L  t  1  (2n+l)n  +  —  I L  —  X I . 


(63)  u(x,t)  =  w  (x,t)  +  2  supf(w  (x*,f)  +  K)  /a  (x')  <  t '  <  t  -  | x  -  X '  )  . 

n  n  n  —  — 

Thanks  to  (61)  and  (62),  u  satisfies  the  boundary  conditions.  Therefore  it  solves  the 
problem  of  the  string  with  an  obstacle  on  Rn  ,  and  the  induction  can  be  pursued.  • 


Ill 


N.rXKK'AI.  SCHEIE 


III. I.  A  numerical  scheme  in  a  backward  cone  for  the  zero  obstacle. 

lot  there  be  given  initial  data  u^  and  u^  on  the  interval  [-T,Tl  .  We  seek  an 

a:  :  mxinat  ion  to  the  :  roblera  p  on  the  backward  cone  T 

,  T  0,T 

t  —  be  a  ate:  ,  and  let  us  define  discretized  initial  data  u!1  and  U*? 

n  0  1 

r»y  t:.e  following  formula 

u''(x>  =  y-  ii  ((|*l)h)  -  u  ([  h)  1  (x-ph)  +  u  (ph)  if  x  (ph,(p+l)h) 


lu-J) 


{;*!)!, 


^  lx)  ^  - 


u  (x')dx'  if 


[ph,  (p+1) h)  . 


TJ.t  correal  ending  free  solution  w'  is  given  by 

ii  i  x+t  i 

v“(x,t)  ■■  -  !  u"  ( x  ♦  t )  1  u’lx-t)  +  r  Ujlx'Jdx*]  . 

x-t 


bet  us  define 
(65) 


h  iU 

2  '  i  2 


for  0 


-h 

Then  w 

i.l 


(66) 


l.j  1  2  :  •  1  2  !  j  -  2  - 

satisfies  the  finite  difference  relation 


•h  ~h  -  h  ~h 

J  -  w.  +w  -  w.  .  ,  . 

i,l  1,1-1  t-1,1  1-1, 1-1 


|i=J 

1  2 


Let  us  define  a  function  u.  on  our  mesh  by 

1,1 

(67) 


*  w'.‘  .  ♦  2  max :  (w'.‘ ,  .,)/»•  i  ,  l '  -  l  ,  i  •  +  i 1 
i,l  i,l  i',l'  -  " 


fi  • . 


We  could  define  j11  al  ternati  vely  by 
1,3 

(68)  u:‘  =  wh  -  2  Kh 

1,1  1,1  1,1 


(6  1) 


•n  -t, 

r  =  min(K  ,  K.  .  ,  ,  -  (w  . )  ) 

l.j  1-1. J  i,1“l  l.l 


if  -n 


botivc  that  K1’  is  not.  the-  disorrti  Nation  of  k  ,  Lot  the*  discretization  of  k.l^, 


w;.t  r»;  I  i  -  t:.i*  :iit  I  ir.  char  act  or  1:- ti  coordinates , 

Iheon  m  ^ .  l.i  t  i.i  m  t:,t  :  ol'ition  of  \j 


-h 


with  zero  obstacle,  and  let  u’.  be 

i .  D 


d  :/•  ( b  7 )  .  Then: 


(7  )  nax  i“  -  i  (  — I  ,  h)  _  '» 

i . ; 

r»  only  on  thr  initial  condition:,. 
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Moreover,  we  have  the  following  bounds  on  the  (approximate)  characteristic  de- 


rivatives : 

(71) 

i  ~h  -h  i 

u.  .  -  u.  _  . 

i  ih 

<  \  |  /  (u 

i/3  1-1/3 

(i-l)h 

i  h  h  i 

.  (-j+l)h 

(72) 

-  2  1  I 

-3h 

1.3  1.3-11 

Proof.  Let  us 

first  evaluate 

w*1  .  -  w(x*  ,■ 

1.3 

,)  (x')dx' 


(u  -  u  ) (x 1 ) d  ‘ 
Ox  1 


.  -  w(x',t')  when  (s',t')  is  in  the  characteristic 
1 . 3 


square  centered  on  (^j^-h  .  ^p-h).  with  sides  of  length  h^2  ,  i.e. 


2l_1  h  <  x'  +  f  <  h  and  ~ -P^1- h  <  x'  -  t'  <  ~2 p--  h  . 

2  —  —  2 


-  2 

w(x',t')  -  Wh  .  =  w(x',t')  -  Wh(  h  ,  iqih) 

1 0  D  z  z 

1 


x.e . 


■  f)  - 

h 

uo 

(ih)  + 

u0(x 

x'+t' 

ih 

/ 

U1 

(y)dy  - 

■  / 

x'-t' 

-jh 

x'+t' 

<-?/ 

'  h 

2 

U  / 

ih 

C 

o 

X 

c 

o 

X 

+  i 

ux)2  + 

(uox 

x'-t' 


tuox  -  ui)2<Jyl1/2i 


“jh 


We  may  then  deduce  from 

(73)  |w(x\t')  -  wh  ,|  <  /%  (  /T  (u2  +  uf)dx)1/2 

1.3  2  jJ,  Ox  1 

that 

|  sup{  (w*?  . )  /i  '  i  .  j '  i  ,  i  + j  0 1  -  sup{  [w(x*  ,t  * ) )  /0  <_  t  •  ■"  t  -  •  x-x*  !  }  |  _ 

^  fh  frT  .  2  2.  .1/2 

<-/  2  (/  (U0x  +  Ul)dx)  * 

-T 

Let  us  notice  that  wh  .  =  w(  ,  — ^-h)#  because  the  approximation  (64)  is 

1,3  2  2 

very  particular. 

This,  in  turn,  gives 
h 


(74) 


(75) 


T  ,  2  2.  .  1/2 

(u  +  u, ) dx) 


max  [  uh  .  -  u  (ijJ-  h ,  h  )  |  <_  iTh  ( /  'u0x  T  “l 
i ,  j  1 ' 3  -T 

This  completes  the  proof  of  (70). 
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We  now  turn  to  prove  (71)  and  (72) . 

_** h  /  **h  ,  i  .  ,  i 

k.  .  »  mintw. ,  .  ,/i'  <  i,  n*  < 
We  can  write  k^  alternatively  as 

(76)  kh  .  =  miniw*!1  .,/i‘  <  i,  j'  < 

i',3  -  ~ 

because  we  know  from  (41)  that  w.  .  >  0 

1.3  - 

the  domain  of  influence  is  not  empty. 
Relation  (76)  implies  that 

(77)  kh  .  =  fn(i)  +  gh ( j ) 

1.3 


Let  us  notice  first  that  if 
j  ,  i'  +  j '  _>  0}. 

j). 

for  i+j  <_  0  ,  as  long  as  we  suppose  that 


where 


~h  _h  . .  h , 
w.  .  =  f  (l)  +  g  (])  , 

1  /  3 

fh(i)  =  min{fh(i*)/i'  £i), 

gh(j)  =  min{gh  ( j  ■ )  /  j  •  £  j}. 

can  be  written  as 

-h  -h  ^  . . .  *h,..,— 

u  =#  +  2|£  (l)  +?  (])1 

i»3  1.3 

not  empty.  If  fNi)  +  gNj)  £  0  ,  then  fh(i-l)  +  gh(j)  £0  ,  and  (71) 
is  immediate.  Suppose  now  that 

(80)  fh(i)  +  gh(j)  <0  . 

We  have  two  cases:  in  the  first  case, 

(81)  fh(i-l)  +  gh(j)  >0  . 

Then,  necessarily 

(82)  fh(i)  =  fh(i)  <  fh(i-l)  <  fh (i-1) 


(78) 
and 

(79) 

Thus,  (68 


if  T  is 


and  thus , 

u?  .  -  uh  .  .  =  fh(i)  +  gh(j)  -  2fh(i)  -  2gh  ( j )  - 
i,3  1-1,3  J 

-  fb "(i-1)  -  gh ( j )  =  -Ifh(i)  +  gh(j)  +  fh(i-l)  +  gh(j)]. 
Thanks  to  (80)  and  (82) ,  we  get 


(83) 


-h 

u.  . 
i,3 


Uj  i  .  |  1  |fh(i)  +  gh(j)|  +  | fh (i-l)  +  gn(j)|  i 

1  if] 

£  IfNi)  -  fh  (i-l)  I  . 
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! 


In  the  second  case. 


if  fh 

(i-l)  =  ?h(i) 

i  -h 

(84) 

u 

1  i. 

-h  „ 

If  f 

(i-l)  >  f  (i) 

fh ( i— 1 )  +  gh ( j)  <  0  , 


Vi  ,  I  1  |fh(i)  -  fh(i-l)  |. 

•L  A#J 


h  **  H  _h  . , .  h  , , ,  ^  ph  #  •  v  A  ~h  ,  , .  .h  . .  _  .  h  ,  , ,  .  ^  h  , ,  _  .  hi  , . 

u.  ,-u.  .  .  =  f  (x)  +  g  (j)  -  2f  (i)  -  2g  (3)  -  f  (i-l)  -  g  (3)  +  2f  (i-l)  +  2g  (3) 
1 *  j  #  D 

=  2fh  (i-1)  -  fh  (i)  -  fh  (i-I)  , 
and,  thanks  to  (82)  we  have 

|uh  .  -  Uh  .  I  <  |fh(i)  -  fh(i-l)|. 

1  1.3  1-1,3 1  “  1 


From  (83),  (84)  and  (85),  we  deduce 

| ir?  .  -  u?  .  .  |  <  |wh  .  -  .  .  |  =  i  |  /  (u  +  u 

'  r,g  i-l ,3  ~  t,3  x-1.3 1  (i-l)  h  0x  1 

The  proof  of  (72)  is  analogous. 


) (x')dx'  . 


We  can  deduce  from  (71)  and  (72)  an  energy  inequality:  let  Iq.Iq  be  9iven  such 
that  -n  <_  i0#i0  n  and  iQ  +  jQ  0  .  Then  we  have: 


(86) 


r  1  I  — h  ~h  I  V  1 

.  ■  h  Ui,j  Ui~l.jn  +  ^  h 

1=-]  +1  Jo  J0  3—1+1 


1  ,-h  ~h  ,2 

u.  .  -  u.  .  .  < 

l0,3  J-o'J'1  “ 


17  f 


-joh 


lu0x|2<x')  +  luj |2 (*' ) Jdx* 


Ill 


periodicity  and  imparity: 


h  ,  .  n 

u  (x)  =  -u  (-x) 
r  r 

for 

-n  h  X 

0  — 

_  0 ;  r  =  0,1; 

u;1(x)  =  -uVx) 
r  r 

for 

nh  •  x  • 

(n«i  )li;  r  -  o,l. 

We  define 

0  #h 

w  by 

f  0  ,h  h 

w  (x , 0)  =  uQ 

(x) 

-noh  - 

x  (n+n^)h 

(88)  | 

,  0,h 

’  ,t  (x,0)  = 

Uj  (X) 

-noh  - 

x  •  (n4nn)h 

r_  0,h 

Cw  =  0 

in  T_ 
n_ 

2 

n  ,  (—  +  n  )  h 

and  let 

(89) 

0  ,h  0,h, 

w.  .  =  w  ( 

1*1 

iiih ' 

¥h>- 

Let 

(90) 

-h  -0,h 

u .  .  =  w .  ,  + 

1*1  1*1 

-0  h 

2  sup) (w.  .  4  K) 

l.J 

/i’  _  i,  j*  _  3,  i'+j*  _ 

where  i 

n+V  j  -  n0'  i+j 

u  . 

Let  us  define  a  subset  p.  'n 

of  2 

•  2  by 

f  i 
(91)  j 

?m'*  =  [n  +  (2m-l)n 

*n  + 

0 

(2m+l) n^] 

'  I-n  *  Om-l)nr  ,  (2m-l)nfJ 

1  (2m-l ) n  ,  n 

0 

+  (2m4l)nr)]  •  [  ( 

!2m-l)n^,  (2ra+l)n,l. 

m.h 

Tnc-  region  P  is  the  discretized  equivalent  (in  i,j  coordinates)  of  the 

region  P  defined  by  (59). 


We  define  'Ja,n  on  the  lower  boundary  of  Rm'^  by 

wm'h  =  un  .  for  i  =  n  +  (2m-l)n  ,  -n  +  (2m-l)n  <  j  <  (2m-l)n 

1,)  1,1  0'  o  -  - 


0  ' 


(92) 


wm'h  =  u*1  for  j  =  (2m-l)n  ,  (2m-l)n  <  i  <  n  +  (2m-l)n 

[  1,3  J  0'  0  -  -  0 

i  wm'h  =  wm'h  for  i  =  (2m-l)n  ,  (2m-l)n  <  j  <  (2m+l)n  , 

i,3  3,i  0'  0  -  J  -  0 


I 


m,h  _  ~m,h 


w. ' .  =  w  ’  .  .  for  (2m  l)n„+n  <  i  <  (2m+l)n  +n,  j  =  -n+(2m-l)n 

1,3  n+3  ,-n+i  —  0  —  —  0 


,  .  m,h 

and  m  R  ,  we  have 

m,h  m,h  m,h  m,h  „  ..  m,h 

(93)  w.  =  w'.+w..,-w.,.,  for  (i,])S(i-l,3~l)  in  R 

1,3  1-1,3  1,3-1  i-l, 3-1 


(94) 


~h 


-m,h 


m  ,h 
R 

m,h 


Then,  we  shall  define  on  Rm'  n  {(i,j)/0  <_n)  by 


m,h. 


u.  .  =  w. ’ .  +  2  sup{ (w  '  . ,  +  K)  /i*  <  i ,  j •  <  j  and  ( i *  , j ' )  c  R  ’  } . 
1,3  1,3  1  ,3  “  “ 


Of  course  (94)  is  the  discretization  of  (63). 

Theorem  9.  Let  u^1  be  defined  by  (93),  and  let  u  be  the  solution  of  on  [ 0 , L] 


with  obstacle  -K.  Then 
(95) 


max  |  u^1  .  -  u(— ^-h,  iU-h)|  <  Cm+'''  /h 

f  j  ^ 

for  (i,j)  in  the  region  R  '  defined  by  (91),  where  C  depends  only  on  the  initial  condi¬ 
tions.  Moreover  we  have  the  following  bounds  on  the  (approximate)  characteristic  derivatives 

i5ifj  •  vi,j  -i  L  L  (uox +  u1><x,>dx,i 


(96) 


(i-l) h 

1  h  h  1  (“j+Dh 

|  ^1,3  ‘  ^j-l1  -21  _{h  (U0X  '  V'x'>d*'l 

if  and  are  extended  to  all  ]R  by  periodicity  and  imparity. 

Proof.  We  shall  replace  the  number  a  defined  in  (87)  by  n  h  ;  for  tnis  new  value  of 


x  ,  we  can  perform  the  construction  of  the  solution  of  as  in  II. 3,  and  we  shall 

m  -m  f  h  ,  ,  m  ,  m  r  h 

compare  w  and  w^  on  the  regions  R  and  R 

Thanks  to  Theorem  8,  the  relation  (95)  is  verified  for  m  =  0  and 
L+  L 


C  '  <2  /  (u^x  +  Ul)dx) 


1/2 


and  the  relation  (96)  is  satisfied  in  R^  . 


Suppose  that  for  a  certain  constant  C  ,  (95)  and  (96)  are  satisfied  in  R 


m-1  ,h 


-2  3- 


(97) 


Then  we  have 

i  ~m,h 
w.  . 
l.D 


-  wm(  i=lh.  i|lh) 


c™  & 


for  i,j  on  the  lower  boundary  of  Rm'h  which  is  the  upper  boundary  of  R™  1,h 


Then  we  have 


(98) 


-  -m,h  m,  i-j  .  i+i  ,  .  i  r  Jn  r 

w.  .  -  w  (  -p-h  ,  h  )  <  5  C  /h 
l,]  2  2  1  — 


m,h 


because  wm,h  ,  (respectively  wm(  ip-h  ,  ii^-h))  is  the  sum  of  at  most  five  terms 
l/j  22 

,  (respectively  wm(  1  h,  1  ^  h ) )  with  i,,j'  on  the  lower  boundary  of 

Rm»h  (respectively  (—  — h,  1  — h)  on  the  lower  boundary  of  R1")  . 

If  we  evaluate  now  the  difference  w”1'^  -  wm(x',t')  when  (x',t')  is  in  the 

1  #  3 

characteristic  square  centered  on  h  ,  1^-  h  with  sides  of  length  hv^  ,  we  have 


(99) 


i  -m.h 
w.  . 

1  1.2 


wm(x,,t')|  1  S(fJh  +  |wm(  i^-h,iiih)  -  wm(x',t') 


but  we  have  for  p  the  equivalent  of  (48),  i.e. 

ljf(x't )l  =  =  l'^'(uox_  V(X_t'0)  i. 

|g(x(t)  |  =  l^(x-t.O)  |  =  | i  (uQx  -  Ul)  (x-t.O)  I. 
if  y^  and  u^  are  extended  to  all  of  ®  by  imparity  and  periodicity. 
Therefore , 


(100) 


.  Ul  . 

w  ~(x,t)  I  =  ^sf  |  (u 


'3C 

liiL 

3n 


Ox  +  V <X  +  t,0)  I  ' 


(x,t) 


1 

72 


(u 


Ox 


u^)  (x-t,0)  | . 


Relation  (100)  allows  us  to  evaluate  w”(^J-h  ,^h)  -  wm(x',t')  , 


2 

L+n 


(101)  I w 


w  (x'tt,)|  <  /2h  (  /  (u  +  u  )dx) 

—  •'Ox  1 


1/2 


Let  us  denote  by  E  the  number 
L+q 


E  =  /  (u  +  u,  )  dx  . 

‘  ox  1 
-a 

Gathering  relations  (97),  (99)  and  (101),  we  obtain-. 

-  u(i^-h  ,  i|ih  )  |  <_  (15Cm  +  2/2E)  Sh  . 
Therefore,  if  we  choose  C=15  +  2v/2E,  we  have 
15Cm  +  2/2E  £  Cm+1  . 

The  proof  of  (96)  is  immediate. 


I  ~h 
u. 

1  ifD 
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Remark  10. 

For 

(i,j)  in  R  '  ,  we  have 

if  >  (2m-l)„o 

and  thus 

1 

+  m  £  ((  i|^-h)/2n0h)  +  j  . 

Therefore , 

if  ( 

,  ~^-h)  converges  to 

(x,t)  as  h  goes  to  zero,  we  have  from 

(95)  : 

t 

2 

1 

'h  ,  i-j  ,  i+j  .  . i  „  2a 

Ui,j  -  U(  2  h*  2  h)>  <  C1 

2  & 

for  all  C 

>  c  , 

,  and  for  all  h  small  enough. 

REGULARITY  IN  SPACES  OF  FUNCTION  OF  LOCALLY  BOUNDED  VARIATION 


This  paragraph  is  dedicated  to  proving  the  following  result  of  regularity  for  an 
arbitrary  concave  obstacle 

1  2 

Theorem  10.  Let  u  and  u,  be  elements  of  H,  (IP)  and  L,  (IP)  respectively, 
-  0  1  loc  loc 


such  that 

du 

(102)  ""dx~  anc*  U1  are  locallV  °f  bounded  variation. 

Suppose  that  u^  and  u^  satisfy  the  compatibility  condition  (22) ,  and  that  the 
obstacle  is  concave. 

Then  for  all  n  ,  the  function 

c,  •••*  —  ( % ,  n ) 

defined  on  l- »", »+«■•)  is  locally  of  bounded  variation,  and  analogously,  for  all  the 


function 


7^  (  V.) 

dri 


defined  on  (-C,+tr)  is  locally  of  bounded  variation. 
Proof.  We  retain  the  notations  of  II. 1: 


U  =  (x/w (x , i (x) )  •  0}  =  u  Ja^,b^[ 


a .  +  t  (a .  ) 
l  l 


b.  +  ;  (b.  ) 
l  i 


-a .  +  ■  (a . ) 
i  i 


-b.  +  ;  (b .  ) 

l  l 


\  Y  <  ; )  =  if 

1 

\  YC)  =  [■!,::.]  if 


i;i-d> 


C  =  (  U  ['.,!))  . 


We  have  the  following  representation  of  the  solution: 


u  ( '  ,  i)  =  f  (  ' )  +  g  (■,)  for 

u(  ,  ■•)  =  f  (  ‘ )  +  g  (-i)  for 


with  the  transmission  conditions: 


fO  *  g(Y(  ))  =  fC)  +  q(Y(  )>  -  ((  -  YC))/  ,1?) 


(105)  f'K)  +  g'(Y(0)  =  -lf’(S)  +  g’(Y(£))l  if  Y  is  one  valued 
and  0  >  Y'  (£)  > 

If  we  differentiate  (104)  with  respect  to  C  on  c  ,  we  get 

(106)  f '(.%)  +  YMUg'lYtU)  =  f'(0  +  Y'  (C)g*  (Y(O)  =  73-* '  <  (C-Y  ( C) ) //2)  (1-Y 1  ( O  ) 
(notice  that  Y  is  decreasing  on  C  ,  and  therefore  almost  everywhere  differentiable). 
For  Y*  (£)  =  0  ,  we  deduce  from  (106)  that 

(107)  f '  (5)  =  f'(£)  =  ;rrv''U  -  Y(O). 

For  0  >  Y’(J)  >  we  deduce  from  (105)  and  (106)  that  f’(S)  +  f  •  (?)  =  *^2  ^  '(S-Y  ( C) 
which  contains  (107).  Therefore,  we  have 

(108)  f1  (?)  +  f  *  (?)  =  /2  v>*  (?  -Y(?)  a.e.  on  C 
and,  differentiating  (104)  on  C  : 

(109)  f •  (?)  =  f '  (?)  a.e.  on  C°  . 

Let  us  denote  by  h  the  function 

(110)  h  ( C)  =  ( C » n) 

where  r,  is  fixed  throughout  the  end  of  this  proof,  and  let 
?0  «  sup{?'  e  X (n)  }  . 

Then,  for  ?  <_  ?Q  , 

h(?)  =  f(?) 


and  ?Q  does  not  belong  to  any  interval  (?^,?p. 


To  evaluate  the  total  variation  of  h  on  a  given  bounded  interval  I  =  [a,b],  we 
have  to  estimate: 


(111)  1 


f  TV(h ;  I )  =  TV  (h ;  I  n  +  TV(h,-C  r>  (?0,+“>))  +  TV(h;  C°  n  (?,+”))  +  I  h  ( +  0)  - 


h(?„  -  0)  1  +  r  (  |  H  { r  .■  +  0)  -  h  ( ? .'  -  0)  I  +  jh(?  .  t  0)  ■  h(C.  -  0)  !)  . 

{i/':0-?i-b;  11  11 


According  to  (108)  and  (109),  we  have: 
f  TV(h;I  (-'-, '0))  +  TV  (h  ;C  (-'Q,  +  ")  )  +  TV(h;CC  r,  (tQ,+»)) 


(112) 


;  ^TV(f;I)  +  TV(/2  v'  >  -  f  ( f ) ;  I ) 

By  hypothesis, 

.reasing,  the  right  hand  side  of  (112)  is  bounded. 


is  positive,  therefore  is  increasing;  as  f  *+  is 
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The  term  |h(q  +  0)  -  h(£Q  -  0)  |  is  bounded,  because  (102)  ensures  that  f  ,  ar.d 

therefore  f  is  locally  bounded. 

The  remaining  term  in  (111)  is  the  sum 

(113)  l  [  |h (C?  +  0)  -  h<£!  -  0)|  +  |h(£.  +  0)  -  h<£  -  0)  |  ]  , 

111 

which  could  possibly  contain  an  infinite  number  of  terms.  Using  (108),  and  (109)  we  can 


write  the  terms  of  (113)  as 

5.-n. 

(114)  |f<q  +  o)  +  f  (q  -  o)  -  Si  *•  (  y-  1  -  o)  |  +  |f(q  +  o)  +  f<q  -  o)  -Jit'  ( +  o)  . 

But  we  have  the  following  inequalities,  deduced  from  the  definition  of  the  line  of  influence 
and  of  the  intervals  [q.CM: 

1  Ci“ni 

f  (q  -  °)  -  -ft  (-^~  -  o)  =  a.  i  o 

1  ^i~ni  + 

f  (£.  +0)  -  4-  *•  (-V^-  +  0)  =  a.  >  0 

(115)  1  ^  1  ~ 

f(q  -  0)  -  (-75-^  -  0)  -  *>"  <  0 

1  ci"ni 

f  <q  +  0)  -  (  yg----  +  0)  =  b*  <  0  . 

We  can  estimate  (114)  by 

[  | a*  +  a" |  +  |b+  +  b‘|  +  Ck(^i+  0)  -  0)1  + 


q  -  n!  q-n’. 

+  [*’  (  •  yg-  -  +0)  -  *•  (  ^  -  0)  1 1 


| a+  +  a. |  +  |b+  +  b. |  <  |a+  +  a” |  +  |b+  -  a+ I  +  |a+  +  a. |  +  'a.  -  b. |  , 

'  1  1  '  '  1  1  ‘  —  '  1  l'  '  1  1  '  1  1  1  1 

and  using  the  sign  conditions  (115) , 

j  |a+  +  a. |  +  |b+  +  b. |  <  2|a+  +  a. |  +  |b+  -  a+ |  +  |b.  +  a. 1  < 

1  1  ‘  ‘  i  l  *  —  ‘  i  l  i  i  i  l  — 


[  <  4  Tv(f(o  -  )»  iq.qi)  . 

Carrying  (117)  and  (116)  into  (113),  we  obtain: 


U/Coiqib) 


l  MS!  +  0)  -  h(S!  -  0)  +  h(c‘..  +  0)  -  h ( f .  -  0)  ,1 
l  i  i  i 


I  4  TV(f»  [i-0,f.5]  )  +  5  TV(  _L  v'  )  :  I-.,.-’']) 


Here,  =  suplf^/S^  <_b}.  The  same  argument  holds  for  the  other  characteristic  de¬ 
rivative.  The  proof  of  Theorem  10  is  complete;  notice  that  we  have  proved,  in  fact, 

that  locally,  TVC^r  (*,n),I)  is  a  bounded  function  of  n  ,  for  all  bounded  1  .  • 

dc. 

Remark  11.  It  is  not  true  that  under  hypothesis  (102),  —■  C * , t )  or  (-,t)  are  of 

bounded  variation  for  all  t  . 

To  see  it,  let  us  consider  the  following  example.  Let 

w(x,t)  =  A-t  -  a(x+t)4  sin  — ~  if  jx  +  t|  <  b 

(119)  •  x+t 

=  A-r  if  lx  +  t|  >  b. 


We  choose  b  such  that  sin  —  =  0  ,  and  a  such  that  the  curve 

(120)  t  =  A -a  (x+t)  4  sin  — 

x+t 

has  always  a  slope  lesser  than  1  ,  for  |  x+t  |  <_  b.  For  this  purpose,  we  differentiate 
(120)  with  respect  to  x  : 

t'  =  4a(l+t')  (x+t) 3  sin  -i—  -  a(x+t)2  cos  — ^r-  *  (1+t')  , 

x+t  x+t 


and  so, 
(121) 


..3  v.2 
4b  +b 

x  2  " 

l-a(4b  +b  ) 


Clearly  |t'|  can  be  made  smaller  than  1  if  a  is  sufficiently  small. 

Then,  we  choose  A  large  enough  to  have 

w (x,0)  =  A  -  ax4  sin  —  >  0  for  |x|  <_  b  . 

duQ  X 

Obviously,  — =  w^(x,0)  and  u^  =  w^tx.O)  are  locally  of  bounded  variation. 

Thanks  to  (121),  the  line  of  influence  is  given  by  (120).  We  shall  now  see  that 

4—  (*,A)  is  not  of  bounded  variation.  The  straight  line  t  =  A  crosses  the  line  of 

n 

influence  infinitely  many  times,  at  the  points 

x  «■  — —  A  for  I— I  <  b,  n  i  TL 
nn  'nr' 

and  we  have 


S  <x’A)  =  -1 

if 

X 

1 

l(2k+2) 

n  A'  (2k+l)i  A!' 

k 

-  0 

or  if 

X 

.  (  1 

1 (2k+l) 

A'(2k+2)  A^' 

k 

-  0 

S  <x,A)  =  +1 

if 

X 

.  (  1 

1  (2k+l ) 

7  _A  '  2k7  "  A)  ' 

k 

-  0 

or  if 

X 

•(set- 

1  -1 

A'  (2k+l)  :t  A)  ’ 

k 

<  0 

This  function  is  not  of  bounded  variation  on  any  interval  containing  zero. 
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V.  CONVERGENCE  of  the  penalty  method. 

Y.l.  Weak  convergence. 

This  paragraph  is  dedicated  to  a  general  (and  unfortunately  coarse!)  study  of 
the  penalized  problem 

f  ...  1  , 

|  -  —  (u,  -v)  -  0 

(122)  J  u  (x,0)  =  u  (x) 

l  '  0 

i  ’u. 

j  — —  (x,0)  =  u  (x) 

where  r  =  sup(-r,0),  and  y  is  an  arbitrary  continuous  function  of  x  ,  and  u^ ,  u^ 
satisfy  the  compatibility  condition  (22). 

Let  us  mention  that  (122)  possesses  always  a  unique  solution;  to  see  is,  it  is 
enough  to  write  (122)  in  the  form  of  an  integral  equation,  and  to  use  Picard  iterations. 

Proposition  12,  We  have  the  following  estimates  for  the  solution  u,  of  (122) : 

b  ^u>  2  3u.  2  :’u,  >u, 

/  [  i^T“  (x,r;(x) )  j  +  |-J£-  (x,o(x) )  \  +  2  (-r^-  )  (x,o  (x)  )  r, '  (x)  ]  dx 


(123) 


-  ^  <l^|2  +  h/jdx 


for  all  Lipschitz  continuous  o  with  Lipschitz  constant  1  such  that 
o(a)  =  o(b)  =  0  . 


(124) 


/  T  <u-.  < 


x',t’)  -  v(x'))  dx'dt*  CfX.t.U^U^) 


x,t 


0  on  (a,b)  , 


where  C  does  not  depend  on  A . 

Proof.  (i)  Estimate  (123). 

We  have  the  identity 

I  1  .  -  3  ’u.  ,u.  j.  lu.,2  ,u .  i2 

I  U‘  U'  v  ^  ^  ,t  3x  ^  it  3x  ^  +  2  "ir  ^  .x  '  +  1  'X  1  + 


-  ( (u  -  y-)‘)2)  =  0  . 


Integrating  on  the  region  i (x,t)/a  ^  x  b  and  0  _  t  _  • (x)  1  ,  we  obtain  the 


identity 


12 


ife?. 

,»■'  •' .  (■ 


,2  ,3uA 


du,  du 

A  A  , 


fb  [| _L  (x,o(x))  j  +  |  — r —  <X,0(x))|  +  2o»  (x)  ^  — T~)  (x,  1  (x)  )  + 

4  t,  uX  '■'t  i’X 

a 

du 

+  i-  ( (u^  (x,o  (x) )  -  v  (x) )  ” )  2  ]  dx  =  /b  (|-+jj-  |  2  +  |u1|2)dx  , 

a 

noticing  that  (u1  (x,0)  -  (x ) )  =  0  for  all  x  .  From  here,  (123)  is  immediate, 

(ii).  Estimate  (124). 

We  integrate  □u>  =  —  (u^  -  ~f)  on  the  backward  cone  Tx  fc: 
x+t  3u 

/  Ou  dx'dt'  =  /  (— -  (x',t  -  |x-x' |)  -  u  (x',0))dx'  - 

A  dt  1. 


*r 

x-t 

Tx,t 

du. 

du. 

o  1 

ft 

<17  (x  +  t  -  f.f)  - 

17  tx  ‘  + 

t'lldt'  = 

3u. 

=  r 

[lt  <X,,t  "  ix  ‘  x'  1) 

+  17  tx’-t  ' 

-  |x  -  x' | ) ]dx'  + 

x-t 

x+t 

3u 

x+t 

+  / 

t  (x'.t  -  |x-x'|) 

'  17  (x*'t  - 

|  x— x ' | ) ] dx '  -  / 

X 

x-t 

Let  o (x' )  =  t 

-  | x-x ' | .  Then 

1 

/  Dux 

(x'  .t'ldx'dt'  <_ 

X+t  dU  dU  X+t 

<  f  [— —  (x',o(x'))  +  o'(x')-T-^  (x*  ,0  (x* ) )  ]dx'  +  /  |u  (x')|dx*  , 

_  3t  3x  „_  +  1 


and  using  Schwarz  inequality  and  (123)  ,  we  obtain 

/  y  (U1  "  v!)  dx'dt'  <_ 

T"  _ 
x,t 

x+t  3u,  3u,  , 


x+t  du,  du  2  ./2  x+r 

/  [— —  (x'  ,0  (x' ) )  +  o'  (x' )— —  (x'  ,o  (x‘ ) ) )  dx’ )  »lr  +  (  /  |  u  (x ' )  |  dx' )  w/2t  <_ 
-  x-t  at  3x  x-t  1 

x+t  ,  du  2  ,  /2 

Y  2v/2t  (  /  |  (  ux  (x' )  | 2  +  (x')|  )dx')  '  .  , 

x-t 

We  need  definitions  of  left  and  right  traces  of  the  characteristic  derivatives  of  a 


function  u  . 


f-i:*  ’ 


’"f 
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The  following  results  were  proved  in  (1):  let  u  be  in  V  (cf.  def .  (17)),  sue:, 

that  Qu  is  a  positive  measure.  Then  the  function 


la,b) 


is  increasing  from  (-a,“)  to  L  (a,b)  for  all  a,b,  and  similarly 

.  ,  du  ■  . 

I  *  —  u,n)  .  .. 

dn  |  [c,d] 

2 

is  increasing  from  I-c,")  to  L  (c,d)  for  all  c,d. 

We  define 

fr  +  h) 


(125) 


~r 

3u 

H 

U,n)  =  lim  - 
htO 

35 

(S,n)=  lim  ■ 

h*0 

dur 

3n 

(n,n)  =  lim 
h*0 

3u* 

dn 

(r;,n)  *  lim 
htO 

dur 

H 

du* 

and  77—  . 

d  S, 

£  « 

dt, 

( 

So 

3u 

3n 


du 

and  for  all  n  larger  than  -£ , ;  analogously,  the  functions  -r—  and 

dn 


dn 


are  defined 


for  all  n  not  belonging  to  the  null  set  and  for  all  £.  larger  than  -r 

Proposition  V.2  and  Corollary  V.4  of  111  tell  us  that 

<•  ,0  (*) )  €  h2.  UR;(1  +  0 '  )dx) 
d£  loc 

(- ,o  (-) )  *  UR  ;  (1  -  o’ldx) 

dn  loc 


du 

3C 


( • ,0  ( • ) ) 


Lloc(lxja(x)  >  Oi,  (1  +  ’ ' )  dx ) 


7-—  ( • , a  ( • ) )  •  L2  ( { x  j  0  (x)  >  Oi,  (1  -  c ' ) dx ) . 

3n  loc 

Notice  that  the  above  traces  are  not  continuous  functions  of  u 
example: 

u  (x,t)  »  1  ♦  —  -  t  if  t<l+i 

n  n  —  n 


We  have  the  followin 


=  t  -  (1  +  -) 
n 


if 


It-. 

n 


ju  1 

Then  -rr-  (x,l)  =  ~rr  ,  V  x,  and 
>  2 


-r  (x,l)  =  - 

i  \  V  z 


Vx,Vn  . 


-32- 


We  may  now  state  the  following  result  of  weak  convergence  of  the  penalization: 


1  2 

Theorem  13.  Given  initial  conditions  u„  <  L,  (IR)  and  u,  '  L,  (ffi  )  such  that 
-  0  loc  1  loc 


u„  >  ^  and  u  >  0  a.e.  on  the  set  {xlu  (x)  =  £(x)},  there  exists  a  function  u 


0  - 
such  that 


1  - 


(126) 

(127) 

(128) 
(129) 


u  c  V 
u  >  •( 

□u  >_  0 

supp  Du  C  {  (x,t)/u(x,t)  =  V?  (X ) } 


(130) 


u  .  l  _  r  ^ 

/  1  |-^r  (x,o  (x) )  |  (1  +  o'  (x) )  +  |-iy—  !x,a  (x) )  |  (1  -  o'  (x) )  )dx  <_ 


3n 


i  !b  <M2 

a 


duo|2 


3x 


)  dx 


i  i 

Jb  il^TT  (*<0  (*) )  |2  (1  +  o‘  (x) )  +  (x,o(x))j2(l  -  j '  (x) )  ]  dx  < 

dt,  d  n  ~~ 


_  du  _ 

I  r  (|u  I2  +  |-^-|2)dx 

b  A 


for  all  Lipschitz  continuous  function  o  ,  with  Lipschitz  constant  1  , 
such  that  j(a)  =  j(b)  =  0  ,  u  •  0  on  (a,b)  . 


(131) 


u (x ,0)  =  uQ(x) 


(132) 


3u 


—  (x,0)  =  u^x)  if  Ug  ( x )  >  v'(x) 


i  3u 


|—  (x,0)  |  <_  u^  (x)  if  Ug(x)  =  »'  (x) 


Proof ■  From  estimates  (123)  and  (124),  we  can  see  that  we  can  extract  a  subsequent 
u  such  that 

(133)  u  -  u  weakly*  in  V  . 

The  weak  *  topology  on  V  is  defined  by  the  semi  norms 


/  u  f^  +  ,/u 


x  r2  1 


'K  f ] 


where  f^.f^  ant*  ^3  are  *n  ( IR+  ; L*  ( IF  ) )  with  compact  support  in  D-  -  I  , 


We  deduce  from  (133)  that 


(134) 


u  •  u  in  C  (IP  ■  Iff  )  with  the  compact  topology. 
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Possibly  with  a  new  extraction 

(135)  —  (u  -  v  )  •  vaguely  in  •  ir  )  the  set  of  measures  on  IR  * 

Therefore 


(136) 


0. 


(  lu 


_  o 

v)  r 


Relation  (123)  gives  a  bound  on 


in  L,  and  thus 
loc 


To  check  (129)  ,  let  (x0<t0'  be  a  Point  such  that  u(x^,t0)  v-  (xf))  ;  thanks  to  (134) 


can  find  a  neighborhood  U  of  (xo,tQ)  and  a  '  0  such  t!lat  u  (x , t )  -  v  (x) 


4  -  *<x0>>  ‘  ‘0  '  V(x't}  '  L’ 

Therefore 


u  j  =  0  for  m 


“0  ' 


and  in  the  limit 


U‘U  =  °  * 


This  proves  (129). 

To  prove  (130),  let  ;  be  given,  and 


be  a  positive  number.  Let  us  define  for 


:  (x)  + 


(137) 


( x )  -  <  x  -  a  +  .  -  . 


,  -x  +  b  -  , 

I  i) 


if 

if 

if 


[a  +  • 
(a  t  • 


o '  b  '  0 


lb  -  0'  b  + 


Then  (123)  implies 


b-.  ..  -u  , 

cl ■  \  \  ,—r  (x,  (x) )  .“  (1+  ■  (x) ) 


—  (x,  ■  (x)  )  "  (1-  •'  (x)  )  ]dx 


b*  •  du 

„ ,  r  0  _ <j  2  2 

)  (  -t —  +  u,  )  dx  . 

'  dx  1 

a- 


Bit  the  left  hand  side  term  of  (137)  can  he  written  as 


b- 


(x)  . 


dx 


(X,t)  (1  <  ’  (x)  )  > 


( x , t )  (1-  ■  (x) )  H 


(X)- 


a.vd  we  can  take  a  weak  limit  in  t:.i  -  d-d  ).  integral  t  i.. 


r  e  (1U). 


-  i4- 


Thus  we  can  rewrite  (138)  without  the  index  u  : 


139) 


Taking 


/  dL  /  0  [||p  (x,o  (x) )  |  2  (l+o' (x) )  +  ll^tx.o  (x) )  |2  (l-o1  (x) )  )dx  < 

-r  '  a+t  o  £  C  C  3n  " 


b+eA  du 

i  If'  +  c">  /  °  !  +  luj  )dx  . 

a‘co 

c '  =  0  in  (139)  and  letting  e"  tend  to  zero,  we  obtain 

b-t  r  r 

/  [  |nrr —  (x,o(x))|  (l+o'(x))  +  (x,o(x))|  (l-o'(x)))dx  < 

..  .  ofj  on 

a+to 


b+t,^ 

dx 


<T°  I2  ♦  |u,|2)dx 


a  -c- 


Letting  go  to  zero,  we  obtain  the  first  relation  of  (130).  If  we  take  c"  =  0 

and  let  e '  ,  then  let  tend  to  zero,  we  obtain  the  second  relation  of  (130). 

The  initial  condition  (131)  is  obviously  satisfied.  It  remains  to  check  (132) . 

For  this  purpose,  let  us  take  in  (137)  o (x)  =0  on  la,b).  Then,  ultimately  we  get 

/b  <x’0)  I2  +  I  ]dx  <_  /b(|-^|2  +  |u1|2)dx. 

a  a 

Using  the  identity 


which  takes  into  account  (131)  ,  we  have 

/a  ||^-(x,0+0)  |  2dx  <  Jb  |  |  2dx  . 

b  a 

As  a  and  b  are  arbitrary,  we  have  eventually 

IIt’  (x, )+) )  I  <  |u  (x)  |  a.e.  on  IR  . 
dt  -L 

When  u^(x)  >  *(x),  we  have  the  first  part  of  (132),  as  locally,  v  -  □  u  =  0  .  ■ 

3uX  3u\ 

We  shall  now  study  the  relation  between  the  strong  convergence  of  ~r —  and  —r— 

dX  d  t 

and  the  verification  of  the  energy  condition  (11) . 

Lemma  14.  Let  be  a  sequence  of  solutions  of  (122) ,  converging  weakly  *  to  a  solu 

tion  u  of  (126)-(132).  Then,  u  satisfies  the  energy  condition  (11)  if  and  only  if 

)Uf  :)U. 

— —  and  — —  converge  to  —  and  —  respectively,  strongly  in  L,  (IR  x  [0,”)). 
yt.  "IX  dt  dX  loc 
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Proof.  Notice  first  that  as 


(u  -  *>) 
y  _ 


•  1  converges  to  v  -1  in  M(1R  *  TP.  ) 
K  K 


weakly,  for  all  compact  K  ,  and  as  (u  -  •f)  converges  to  zero  uniformly  on  compact 


sets ,  then 
(140) 

for  any  compact  set  K 


/  [((  u  -  vO ”) "Vy ) dx 1 dt 1  -*•  0 

'  JJ 

K  ^ 


Let 

a. 

(x'l 

1  =  h 

-  |  X— X 1 

h 

x+h 

r 

r  1  3v  , 

F 

dh 

[  —  (x,< 

0 

x-h 

1 3?  ’ 

(141) 

t+x-x' 

= 

F 

dx' 

/  2 

x-t 

0 

If  the  limit 

of 

the 

sequence 

x+h 

du„  2 

(142) 

/ 

dh 

/ 

<1-5-1 

0 

x-h 

The  value 

of 

(141)  for  v  =  i 

x+h 

dUQ  2 

(143) 

F 

dh 

I 

( I'dx'l 

0 

x-h 

3C 


x+t  t-x+x’ 

/  dx'  / 
x  0 


||^<x.,f)|2dt. 


0<t  *  <t-  x-x* 


((u  -  *f)”)2dx'dt 

y 

3u 


3U 

Tn 


And  according  to  (140),  the  limit  of  (143)  is  (142).  Therefore,  as  (resp.  -7^7  ) 

converges  weakly  to  -j^r  (resp  in  L^^UO,")  x  ®+ )  >  and 


3u  ,  3u  _  ,  .  ...  , 

lim  ^  lTf'|2dx'dt'  +  ^  l“arl  dx'dt'  =  J  I  if  I  dx'dt'  +  /  |j-|  dx ’ dt 1 

M+0  A  .  H  B  „  3,1  A  .  H  ^ 

x,t  x,t  x,t  x,t 


where  A  ^  -  {(x'.t')  e  T  ./x*  <  0}  ,  B  ,  -  T  .  \  A  .  r  we  can  conclude  that  the 
x,t  x,t  —  x,t  x,t  x,t 

9u  3u  .  . 

ir  P  J  P  9u  _  dU  .  . 

convergence  of  — r  and  — —  to  —  and  —  is  strong. 

at,  9n  oc,  on 

3u  3u 

Conversely,  if  -r~  (resp.  — —  )  converges  strongly  to  —  (resp.  — ) ,  then  it  is 

at,  on  o  c,  v*  n 

straightforward  to  pass  to  the  limit  in  (11).  • 
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This  problem  can  be  solved  immediately  and  qiv.  .  .  . 

sequence  u  is  u  ,  and  lemma  1-;  allows  as  to  I  . 

Remark  16  .  Suppose  we  replace  the  function  r  :  y  .1 


is  continuous,  strictly  lien  a:-  c.  ;  or,  i-  ,  J 


then  the  penalised  j  rot-lem 


- .  (u  -  >') 


u  (x,  1)  -  u  (x) 


— —  (x,o)  =  u  (x) 


can  oe  studied  as  above;  we  get  Theorem  14  with,  •ilr.vr.se  no  change  in  the  proof.  ",oi<  c  r 

a  phase  plane  analysis  shows  easily  that  in  the  case  of  initial  data  (144)  tne  limit,  of 

u  is  the  function  (146).  vie  chose  the  specific  penalization  (122!  because  of  its  sim;  lic- 
ity.  We  need  an  integral  solution  of  the  linear  Klein-Gordon  equation  with  initial  values 
given  on  a  curve  t  =  : (x)  .  It  is  the  object  of  the  next  lemma. 

Lemma  1 7 .  Let  w  be  a  solution  of  the  wave  equation  on  the  set  S  =  ■  (x,t)/:(x)  _  t  ■_ 

t  -  j x— Xq j  >  where  is  a  Lipschitz  continuous  function  with  Lipschitz  constant  1  . 

Then  the  unique  solution  on  S  of  the  p.roblcm 

‘  1 

_.u  *  —  u  =  0 

i 

(147)  I  u(x,-(x))  -  w (x, r (x) ) 


—  (x,  •  (x) )  =  —  (x,  (x))  a.e.  on  •x/l',(x)1  1 


is  given  by 


u  ( x ,  t )  =  w  ( x ,  t )  - 


r  J,  (’ 


(t-t1)  -  (x-x 1  )  ‘ 


■)w(x*  ,  t 1  )dx'dt'  , 


is  the  Bessel  function  J  . 


Proof.  We  verify  that  if  w  is  a  solution  of  the  wave  equation  in  the  whole  plane 


and  if, 


then , 


w(x,t)  = 


0  if  t  <  o(x) 

w(x,t)  if  t  >  o(x) 


lCw,v5>  =  -  /w  (x,a  (x) )  (^t  (x,o  (x) )  +  a*  (x)^x  (x,o  (x) ) )  dx  + 
+  /[  wfc(x,  (x) )  +  o’  (x)wx(x,a(x) )  ]^(x,a(x))dx. 


Solving  (147)  amounts  to  find  a  solution  of 

1 ... 

|  {  !x,t)/t  >  o  (x)  } 


(Du 


i) 


U(X,t> |{ (x,t)/t  ^  0(x) }  =  0 
u(x,o(x))  =  w(x,a(x)) 

(x,a(x))  =  (x,o(x))  a.e.  on  {x/|o'  (x)  |  <  1} 

ofc  ot 

which  can  be  written  as 
(ISO)  u  =  w  -  y-  fi*u 

where  8  is  the  elementary  solution  of  the  wave  equation  defined  by 

if  t  _>  x 
j^O  elsewhere 

The  convolution  equation  (150)  has  a  unique  solution  given  by 


£(x,t)  = 


(151) 


l 

n=0 


(-1) 


( 


8)  *  w  . 


By  a  simple  inductive  calculation  in  characteristic  coordinates,  we  obtain: 


(152) 

Therefore 

(153) 


(  *  6  )  (€,n)  =  6  )n_1  - - - T  • 

k=l  ( (n-1)  ! ) 


<*>  n  n  /  2 

l  ~r<k: i8  x^)  •  J0(/i-T- 

n=l  A 


Together  with  (153),  formula  (151)  gives  (149).  ■ 

We  can  now  state  the  theorem  for  convergence  for  penalized  solutions: 

duQ 

Theorem  18.  Let  uQ  and  u1  be  such  that  — and  u^  are  locally  of  bounded  variation 
(154)  and  suppose  they  satisfy  the  compability  condition  (22).  Then  the  solution 

u  of  (122  converges  to  the  solution  of  when  >  goes  to  zero. 
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proof.  Let  us  first  notice  that  on  I  ,  the  complement  of  the  domain  of  influence,  we 
have,  if  u  is  the  solution  of  P  : 

[  Du  =  0  ] 


Therefore 


u^  =  u  on  I 


and,  in  particular, 

u^(x,t(x))  =  u(x,t(x))  =  w(x,t(x)) 

(155)  ■  — i  (r ,  t  (x) )  =  1^-  (x,  t  (x)  -  0)  »  Ijr  (x ,  t  (x) )  a.e.  on  {x/|t'(x)|  <  1; 

3r  or 

where  we  recall  that  T  ,  the  line  of  influence,  is  Lipschitz  continuous,  with  Lipschitz 
constant  1  . 

We  shall  now  use  assumption  (154)  to  obtain  more  information  about  the  line  of 


influence.  We  need  the  following  notations: 


5^  =  {  (x,t)/x  2L  1 1 1  } 
Q2  =  {  (x,t)/t  >_  |x|  ) 
Q3  =  {(x,t)/x  <_  -  1 1 1  } 
g4  =  {(x, t)/t  _<  —  l X I  > 


We  shall  denote 


<x,t;Q.) 


lim  —  (x+h,t+k) 
(h,k)-0  ' 

(h , k)  •  Q. 


3w 


Thanks  to  (154)  ,  —  (x,t;Q.)  is  defined  for  1  i  <  4,  and  we  have  the  formula 

dt  i 

(158) 

with  notations  (125).  We  have  analogous  formulae  for  the  three  other  limits. 


If  -V?  [JSf  +17T  (x'tn  ' 


Lemma  19.  Let  x  be  such  that  t*(x)  is  defined  and  |t'(x)|  <  1  .  Suppose  that 
9w 


(159) 


max  t r*  (x,t;Q.)  <  0  . 

,  .  dt  1 

l<i<4 


Then  there  exists  a  neighborhood  (x-e,x+e)  of  x  such  that  |x'-x|  <  e  =*  i'(x')  has 

left  and  right  limits  at  every  point  and  |t'(x'  ±  0)  |  <  1  ; 

sup  )  <_£<() 

l<i<4  3t  '  ,TlX  >,Qi'  -  *  U  • 

Proof .  The  hypothesis  (159)  implies  that,  in  a  neighborhood  N  of  (x,t(x)): 

SUP  (x.  t.,0  ,  <  <  o 

l<i<4  3t  lX  ,C  'V  -  *  0  ' 

therefore  w(x‘  ,«)  is  strictly  decreasing  for  x’  close  enough  to  x  ,  and  moreover 

if  k  is  so  chosen  that  w  (x,x(x))  +  k  w  (x,t(x))  <  0  ,  then 

t  x 

w (x  +  kh , t  (x)  +  h)  <  0  . 

Thus,  there  exists  a  unique  solution  to  the  problem: 
w(x’  ,a'  (x1 ) )  =  0 

(16°) 

max ( |x-x' | , | a (x' )  -  t (x) | )  ^  a  where  a  is  a  small  positive  number. 

To  prove  that  o  is  identical  to  T  in  an  interval  [x-a’,x+a']  where  a'  may  be 
smaller  than  a  ,  we  have  to  check  that 

|c'(x')|  <  1  a.e.  on  [x-a',x+a’]  . 

The  function  a  is  continuous  indeed,  as  w  is  continuous  and  t'  =  o(x')  is  the 

unique  solution  of  w(x',t')  =  0  in  N  .  We  may  not  directly  differentiate  the  rela¬ 

tion  w(x',o(x)'))  =  0  ,  as  we  do  not  have  the  assumptions  of  the  implicit  function 

theorem.  But,  with  the  very  same  argument  as  in  this  theorem,  and  using  notation  (157) , 
3u 


and  its  analogue  for  q—  ,  we  have 

oX 


(161) 


w  (x'+h  ,a  (x+h) )  =  w(x,o(x'))  +  w^  (x '  ,o  (x' )  ?Q  )  h  + 

+  wt(x',o(x');Q2)  (a(x'+h)  -  o(x'))  +  (  |  h  |  +  |o(x'+h)  -  o(x')|) 

for  all  h  such  that  (h,o(x'+h)  -  o(x')>  . 
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Here 


is  a  function  such  that 


lim 

k-0 

By  a  standard  argument 
(162) 


<k) 


lim 

h-0 

(h  , o  (x'+h)  (x1 ) )  ■  Q. 


'  o(x' 


+h)  -  ,-(x') 


w  ( X '  ,  . ( X • ) ;Q, 
x  1 


J  W  (x'.-(x' 


The  same  result  holds  in  the  three  other  quadrants  Q^,  and  and  hy  choosing  1 

adequately  small  we  shall  have 

j  wx(x',o(x')»Ci)  j 


w  <x\j<x*);Q.)  1-e 

t  — 


for 


|x-x' ! 


and  thus 


(h,j  (x'+h)  -  i(x'))  ■  for  h  >  o  ,  small  enough 

(h,  (x'+h)  -  (x  * ) )  ■  Q3  for  h  <-  0,  jh|  small  enough 
i  a  •  (x  * )  j  ^  1-  ■  a.e.  on  (x-i,,x+<'l 
t(x')  =  o(x'),  and  has  right  and  left  limits  at  all  points  of 

(x~:*',x+u'J.  • 

Let  us  compare  locally  the  solution  of  the  linear  Kloin-Gordon  equation  (147)  to  the 
solution  of  an  approaching  problem  with  simpler  initial  data.  Lot 


iUQ) 


0  ’ 


1o(xo)  =  m 


WV 

wt(xo’to)  ' 


T0(x)  =  *0  +  m(x-xo)( 

„0(x.t)  =  wt(x0,t0)(t-t0)  +  WX^X0'^0^  ^X-X0^  ' 
u0(x,r0(x))  =  w0(x,t0(x)>  =  0 

JU  -W 

—  (X, :„<*)>  -  —  (x,i0U))  =  wt(x0,t0) 

SQ  =  { (x,t)/t  : 0  (x )  J. 


Then : 


t-Vro(x*xo) 


u,  (x,t)  =  /  (1-m  )  w  (x  ,t  )  sin  ■  .  '  y 

0  fc  °  0  /  (1-m2) 
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With  the  help  of  (148),  we  have 


(163) 


u(x,t)  -  uQ(x,t)  =  w(x,t)  -  WQ(x,t) 


.  _L  r  J  f 

2A  T-J  0^ 
x,t 


4t-t‘  )2  -  (x-x1)2 


1  (w  •  1  )(x*,t*)  -  (w  •  1  )  (x*  ,t')  Jdx’dt' 

b  U  bQ 


ix-x0|  +  |t-t0|  <  C  A 


Let  us  estimate  (163)  for  x  and  t  such  that 

(164) 

and  under  the  hypothesis  that  |i*(x  )|  <  1  and  that  w^tx^t^)  anc*  Wx*X0,t0^  are 
well-defined.  Then: 


|w(x,t)  -  wQ(x,t)|  <^o(|x-xQ|  +  |t-tQ|)  =  o(^A). 


To  estimate  the  integral,  let  us  first  note  that 


w  •  1  -  w„  •  1  <  w  -  w  I  •  1 

1  S  0  Sq 1  —  1  0  S  n  SQ 


This  relation  comes  from  the  fact  that,  locally,  w  •  1  =  -w  and  w^ 


1  =  -w 

S0  ° 


We  define  new  variables  X  and  T  by 

t-t1  =  t/a 

X-X*  =  x/i  . 


Then  the  integral  expression  in  (163)  is  estimated  by: 

A  /  |  J  </t2-X2)  |  |w(x-x/A,t-T/A)  -  w  (x-x/I,t-T/i)  |l  (x-X/A  ,t-TfA) dXdT. 

T+  °  0  Is  "  s0 

0,0 

But: 

|w(x-x/A  ,t-T/X)  -  w0(x-x/A,t-T^)|  o  (  I  x-x/a-Xq  |  +  |(t-T/A-tQ|) 
and  we  have  to  check  that  {(X,T)  <  T*  /  (x-xA ,t-T/A)  <  S  "  is  bounded.  This  set 

can  be  written  as : 

f  j-T(x-x/A)fT  (x-x/I)  \\ 

\  (X,T)/|x!  T  <  —  -  min  ]  -  \ 

l  i/i  (  /\  i  i 

and,  using  the  fact  that  |  t  '  (x  )  |  <  1,  this  set  is  bounded,  under  the  condition  (164). 

Thus,  immediately: 

(165)  |u(x,t)  -  uQ(x,t)|  =  otA). 

A  consequence  of  (165)  is  that,  for  >  sufficiently  small,  the  solution  u  of  (147)  is 
negative  on  the  set 

(166)  T  , - —  {  (x,t)/t  •_  :(x)!. 

*0,t0  +  )/•  > 
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This  uses  the  fact  that  ufc  <  0  on  a  neighborhood  of  xQ  ,  as  was  proved  in  Lemma  15. 

Therefore,  on  the  set  (166) ,  the  solution  of  the  penalized  problem  (122)  is  the  solu 
tion  of  the  linear  problem  (147),  for  X  small  enough.  We  have  thus,  for  (x,t)  on  th 


set  (166)  : 

3u 


—  (x*  t'  )  =  —  (x 

3t  1  ,c  1  at  ' 


•.f)  -  I  H  (x*,t- 

(x* ,t-]x-x' | ) eS  3 


x-x'  )dx' 


1  r  T.f  /(t-f)2  -  (x-x1)2  ) 

xxr’  4 - H  J 

x,t 


)2  -  (x-x'-2 


A 


(t-t')(w  •  ls> (x,t)dx'dt. 


Reasoning  as  for  (165),  we  can  prove,  under  assumption  (164)  that 


3ux  3uo 


at  at 


<X,t) 


=  0(1)  , 


au 


at 


X  t-t  -m(x-x  ) 

-  (x,t)  -  w  (x  ,t  )  cos 


t'  O'  0 


and,  in  particular 


/x  (1-m2 


=  0(1) 


(167) 


lim  3uX  (x  ,t  +  (H-e) /x  (1-m2)  =  +w  (x.,t.)  cos(ir-e) 

X-0  at  0  0  too 


on  ( Xq/t 1 (xq )  <  1  and  wfc (xq,t (xq) ;Q. )  <  0  ,  i  =  1,...,4}. 


Analogously 

3ux  / - ~ 

lim  -r—  (x  ,t  +  (rr-c)/X(l-m  ))  =  w  (x„,tn)  cos(tt-e) 
X+o 


(168) 


X  0’  0 


on  (Xq/t*  (xq)  <  1  and  wt(x0,T (xQ) s >  <  0  ,  i  =  1,...,4}, 


and  (167)  and  (168)  in  turn  imply: 

3ux  / - - 

lim  (x0,tQ  +  (n-e)A(l-m  ))  =  w  (x0,tQ)  cos<lr“e>  , 

A-*0 


(169) 


1  3«a  ) — r~ 

Um  T7T  u0't0  +  (,,”e)  A (l-m  ))  =  w_(x„,t„) 
X-K) 


n  o'  o 


cos  (h-e)  . 


Therefore,  the  limit  u  of  u  satisfies: 


(x,t  (x) )  =  -  (x,t  (x) ) 


a.e.  on  {x/| t • (x) J  <  1). 

This  proves  that  u  is  indeed  the  solution  of  P  . 
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